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1. Introduction 

When a mesoscopic wire made of normal metal N is in contact with a su- 
perconducting reservoir S, Andreev reflection (AR) occurs [1]. This affects 
the electronic properties of the wire [2]. In this article we address both 
experimentally and theoretically the following question: how does Andreev 
reflection manifest itself in shot noise measurements, and what physics can 
we deduce from such measurements ? Our discussion will rely on high fre- 
quency measurements performed on various NS structures, that can be 
found in refs. [3, 4, 5]. 

The shot noise is a direct consequence of the granularity of the electric 
charge. Even though an electric current exists only if there are free charged 
carriers, the value of this charge does not directly affect how much current 
flows through a sample when biased at a finite voltage. One has to investi- 
gate the fluctuations of the current to determine the charge of the carriers. 
As a consequence the measurement of the shot noise offers direct access 
to the elementary excitations of any system, through the determination of 
their effective charge. An effective charge different from le directly reflects 
how, due to their interactions, the electrons are correlated, as in a super- 
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conductor, a 2D electron gas in the fractionnal quantum Hall regime or a 
ID Luttinger liquid. 

In an NS system, electrons can enter the superconductor only in pairs. 
Thus, the elementary charge participating in the electric current is no longer 
e, but 2e. Thus, a naive expectation is that the shot noise should be doubled 
in the presence of an NS interface, as compared to the case of a normal 
metal. In the following we explore this expectation, to see how Andreev 
reflection affects current noise. In particular we show how a deviation of 
the effective charge from 2e reflects the existence of correlations among 
pairs of electrons. 

In equilibrium the current noise power Si is given by the Johnson- 
Nyquist formula, which relates the current fluctuations to the conductance 
of the sample G: Sj(V = 0,T) = 4k B TG, where T is the electron temper- 
ature. Whether the fluctuating current is made of single electrons or pairs 
affects equilibrium noise only through the conductance, which may depend 
on the state (N or S) of the metal. This picture is valid at low frequency 
(hco <C k B T). It is the 'classical' regime. At finite frequency uj a quan- 
tum mechanical treatment of noise is necessary. The general expression for 
the equilibrium current noise is given in terms of the reduced frequency 
w = hw/{2k B T) by [6]: 

Sj{V = 0, T, uj) = Ak B TG(T, uj)g(w) (I) 

where the function g(x) is defined as: g(x) = x coth x. As in the classical 
regime, all the physics is contained in the conductance G, but here G is 
the real part of the complex and frequency dependent admittance of the 
sample. The g(w) term accounts for statistical distribution of an excitation 
of energy fiuj at thermal equilibrium. Through g(w) the finite frequency 
adds an energy scale at which the classical-to-quantum noise crossover takes 
place: fiuj = 2k B T, or w = 1. 

Since we are interested here in shot noise, let us discuss first the case 
of the normal tunnel junction. The shot noise of the tunnel junction at low 
frequency is given by [7]: 

Sj unnel (V,uj = 0, T) = 4k B TGg(v) (2) 

with the reduced voltage v = qV/(2k B T). Here g(v) interpolates between 
Johnson noise (g(0) = 1) and shot noise (g(v 3> 1) = v). Hence the charge 
q appears at two levels: it can be measured through the equilibrium-to-shot 
noise crossover occuring at qV = 2k B T, or through the magnitude of the 
noise at high voltage, such that Sj unnel = 2ql. At finite frequency, the noise 
emitted by a tunnel junction is: 



Sj unnel (V,io,T) = 2k B TG{g{v + w) + g(v - w)) (3) 
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Thus, a finite frequency investigation offers another way to measure q, 
through the classical-to-quantum noise crossover occuring at qV = huj. 

In the case of a diffusive mesoscopic wire, the shot noise is reduced 
by the disorder, but the principles above are still valid. Thus, AR should 
show up in the shape of the noise spectrum as a doubling of the classical- 
to-quantum noise crossover frequency, occuring at huj = 2eV. To perform 
such a measurement as a function of frequency, one needs to have a precise 
knowledge of G(u>) and of the frequency response of the experimental setup 
at high frequencies. (For example, T = lOOmK corresponds to V = 8.6/uV 
and w/27r = 2.1GHz.) As a consequence, the direct measurement of the 
frequency dependence of the noise spectrum at fixed voltage has never 
been accomplished. The noise measured in a narrow frequency band is 
expected to show the same crossover, but as a function of the applied 
voltage. This is a much easier (but still difficult) experiment, which we 
shall report in section 3. The measurement has been performed so far on 
a normal metal wire, but could also be carried out on an NS sample. In 
Section 4 we report measurements of photon assisted noise in an NS wire, 
which provide an alternative to the measurement of the crossover frequency. 
In this experiment a high frequency excitation is applied to the sample. The 
shot noise develops features as a function of voltage V each time qV is a 
multiple of the energy of the incident photons hu>. For the NS wire, q = 2e. 

The discussion above treats the consequences of the AR due to the 
energy qV, as compared to fegT or Two. It investigates the effect of AR 
on the distribution statistics of the electronic excitations (which involve 
pairs of electrons) rather than the effective charge they carry. Specifically , 
the phenomena that have been measured and discussed above are related to 
steps in the distribution function, as will be discussed in section 6. A better 
measurement of the effective charge is in the fully developed shot noise 
regime (eV 3> hsT). Here the noise is determined by the fact that electrons 
are paired and also by the interferences and the correlations that can exist 
among pairs. In section 5 we report measurements of phase dependent shot 
noise in an Andreev interferometer, which point out such a sensitivity of the 
effective charge to pair correlations. Section 6 is devoted to the theoretical 
investigation of the effective charge deduced from the shot noise. Section 
2 contains information about sample preparation and experimental setup, 
common to the experiments reported in subsequent sections. 

2. Experimental considerations 

The experiments we report in the next sections use samples prepared with 
similar methods, and measured with similar detection schemes. Each sam- 
ple consists of a metallic wire or loop between two metallic reservoirs, either 
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normal or superconducting. The measurements are performed through con- 
tacts to the two reservoirs. This allows dc characterization of the sample 
(two-contact differential resistance Rdiff = dV/dl, measured at ~ 200Hz) 
as well as high frequency measurements. Even for the measurement of the 
effective charge of the Andreev interferometer, which does not intrinsically 
call for the use of rf techniques, high frequency measurements have been 
chosen for their extremely high sensitivity (the signal-to-noise ratio is pro- 
portionnal to the square root of the bandwidth times the integration time 
of each measurement). 

The measurements were performed in a dilution refrigerator at a mixing 
chamber temperature T ~ 50 mK. At low temperature the electron energy 
relaxation is dominated by electron-electron interactions [8] and the asso- 
ciated inelastic length L ee is larger than L, so the transport in the device 
is elastic. We have not conducted weak localization measurements on these 
samples. From these, the phase coherence length could have been ex- 
tracted (L ee and L v coincide if the dominant phase relaxation mechanism 
is electron-electron interaction, which is likely in our samples, otherwise 
L ee > Ltp) [9]. Nevertheless, we observe significant harmonic content of the 
R vs. flux curve of the interferometer (data not shown). The n th harmon- 
ics decays as exp(— nL/L*(T)), where L is the distance between the two 
reservoirs. The empirical characteristic length L* includes phase beaking 
mechanisms (L^,) as well as thermal averaging (usually described by the 
thermal length L T = (HD /k B T) 1 / 2 ). We obtain L*{T = 50mK) = 800nm 
of the order of Lt, which implies that 3> Lt, i.e., L v is much larger 
than the sample size L = 540nm. Other measurements on a longer inter- 
ferometer (L ~ 1/xm) also gave evidence that L ee > L in that device. The 
wires described in sections 3 and 4 are shorter and thus are also likely in 
the regime L < L^. 

2.1. SAMPLE PREPARATION 

The samples studied have been patterned by e-beam lithography. All are 
made of thin (lOnm) evaporated gold wires between thick metallic reser- 
voirs. The N wire and N reservoirs are deposited using a double angle 
evaporation technique [10] in a single vacuum pump down. Sputtered Nb 
is used for the thick (80 nm) S reservoirs. The transparency of the NS in- 
terface has been achieved by ion beam cleaning before Nb deposition. We 
estimate that the interface resistance is less than 1/10 of the wire resis- 
tance. Theoretical calculations which consider this extra resistance show 
that its effect is negligible. The Au wires have a temperature independent 
sheet resistance of the order of ~ 10 Q per square. The Au reservoirs are 
70 nm thick and have a sheet resistance of less than ~ 0.5 Q per square. 



Figure 1. Schematics of the different samples that have been measured, (a) N wire 
between N reservoirs (L = 200nm, D = 40 cm 2 /s). (b) N wire between N and S reservoirs 
(L — 280nm, D = 30cm 2 /s). (c) Andreev interferometer (L = 540nm, D — 33cm 2 /s). 
(d) SEM picture of the Andreev interferometer; N contact on right not shown. The scale 
bar corresponds to 1/im. 



2.2. EXPERIMENTAL SETUP FOR HIGH FREQUENCY NOISE 
MEASUREMENTS 

The experimental setup we used to perform the noise measurements on the 
Andreev interferometer is depicted in fig. 2. The current fluctuations Si in 
the sample are measured in a frequency band Af from 1.25 to 1.75GHz 
using an impedance matched cryogenic HEMT amplifier. The noise emit- 
ted by the sample passes through a cold circulator, employed to isolate the 
sample from amplifier emissions. It is then amplified by the cryogenic am- 
plifier and rectified at room temperature after further amplification. The 
detected power is thus given by Pdet = GAf (kBT out + UbTa) where G is the 
gain of the amplification chain, Ta ~ 6.5K is the noise temperature of the 
amplifier, and T ou t is the effective temperature corresponding to the noise 
power coming from the sample {T ou t = 0.04 — 0.6K for V = — 150 //V 
for the case of an NS interface). We determine GAf and Ta by measuring 
the sample's Johnson noise vs. temperature at V = and its shot noise at 
eV » (ksT, Ec). (E c = HD/L 2 is the Thouless energy of a diffusive wire of 
length L; it is the energy corresponding to the inverse of the diffusion time 
along the wire). We modulate the current through the sample to suppress 
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Figure 2. Experimental setup used for high frequency measurements. The inner dotted 
line correspond to the mixing chamber at T = 50mK, the outer one to the He bath or 
the 4K stage of the dilution refrigerator. The position of the coupler and the circulator 
are those used for the measurement performed on the Andreev interferometer only. 



the contribution of Ta- We measure dPdet/dl. This gives dT out /dI. T out is 
related to the noise emitted by the sample through: 

T ou t = (1 - \r\ 2 )T N + \T\ 2 T in (4) 

and Sj is given by Sj = AksT^ReZ^j:. Here T/v is the sample's noise 
temperature, Z^ff is the complex differential impedance of the sample at 
the measurement frequency, T is the amplitude reflection coefficient of the 
sample and Tj n the external noise incoming to the sample. In eq. (4), the 
first term on the right represents the noise emitted by the sample which is 
coupled to the amplifier. The second term represents the external noise the 
sample reflects. 

In order to determine Sj at finite frequency, it is necessary to know 
both Zdiff and T at the measurement frequency. Z^iff is deduced from the 
measurement of T through: 

1 ' Z diff (u) + Z (u) { ) 

where Zq the impedance of the measurement apparatus. Zq is ideally real 
and equal to 500. In practice it has an imaginary part and is frequency 
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dependent (due to finite return loss of the amplifier or isolator, parasitic 
capacitance in parallel with the sample, inductance of the wire bond, etc.). 
Thus a careful calibration is necessary to have a reliable measurement of 
Zdiff[H]. However, since our samples have a resistance close to 50f2, the 
amplitude of |T| 2 is of the order of a few percent, and can be neglected in 
the first term of eq. (4). This is not always the case for the second term. In 
the measurement performed on the normal wire (section 3), no circulator 
Wcis used . T^yi 

30K » T N for this broadband (20 GHz) amplifier. The 
impedance of the sample (a very short gold wire) is voltage and frequency 
independent, so that the noise reflected by the sample adds up to the total as 
a voltage independent (but frequency dependent) constant. (Tj n depends on 
frequency because the amplifier emission does). For the NS wire (section 4), 
a circulator placed in liquid helium has been used. The circulator attenuates 
by 20dB the noise emitted by the narrowband amplifier towards the sample 
(Temit ~ 2K). In that case Tj n is equal to the temperature T = 4K of the 
50f2 termination of the circulator. Tj n = (T em j t /100) + AK > T em i t , but 
Temit drifts, which adds drift to the measurement. For this experiment one 
was interested only in the voltage dependence of the features of Si, the 
corrections due to Tj n were not significant. 

For the precise measurement of the effective charge in the Andreev in- 
terferometer, the magnitude of Si is of interest; it is thus crucial for T; Ln 
to be minimized and stable. For this experiment we therefore placed the 
circulator at T = 50mK. We also measured the variations of |T| 2 by send- 
ing white noise to the sample through the unidirectionnal coupler (see fig. 
2) and detecting the change in the noise power. This determines the rela- 
tive variation of the amplitude of the reflection coefficient (as a function of 
voltage or magnetic flux) over the bandwidth used for the noise measure- 
ment. We draw the following conclusions: i) the variations of |T| 2 are small 
enough to be neglected, allowing us to take T = in the data analysis. 
This is confirmed by the fact that at V = 0, where T/v = T, Rdiff and T 2 
are flux dependent; yet T out at V = 0, does not depend on the flux (see 
eq.(4)); ii) the impedance of the sample at the measurement frequency is 
different from its dc value. This conclusion is seen from the following : if 
Zdiff(uj) were equal to its dc value Rdiff = dV/dl, then whatever Zq{u) 
is (i.e. whatever the imperfections of the experiment are), |T| 2 plotted as 
a function of Rdiff should collapse into a single curve for all the values of 
flux and voltage. As shown on fig. 3, this doesn't occur. This means that 
the impedance of the sample is not equal to Rdiff = dV/dl measured at 
low frequency. It has an flux- or voltage-dependent imaginary part, or its 
real part is not simply proportionnal to Rdiff - This observation deserves 
more study, through measurements of the amplitude and phase of the re- 
flection coefficient as a function of flux, voltage and frequency. However, 
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Figure 3. Measurement of the variations of the power reflection coefficient |r| 2 vs. dc 
resistance of the Andreev interferometer. Each curve corresponds to a fixed magnetic 
flux and a varying current, between — 2^A and 2/xA. The reflection coefficient has been 
arbitrarily rescaled to its value at V = 0,3> = 0. Conclusions are given in the text. 



this effect is small, and for our present study of noise, we simply use Rdiff 
to determine Si from T/v • This is also justified by the fact that in our short 
phase-coherent samples transport is elastic. Thus, the ac conductance is 
given by the dc /(V) characteristics shifted by ±fruj/e ~ ±6/uV [7]. Since 
the characteristic scale for changes of dV/dl is ~ 30/xV, finite frequency 
corrections to Rdiff should be small. 

In the experiment described in section 3, the noise needs to be measured 
over a broad frequency range. Thus, a broadband (1 — 20 GHz) cryogenic 
amplifier has been used, even though it has a higher noise temperature 
(Ta ~ 100K) than a narrow band amplifier. Also a circulator cannot be 
used. The noise at different frequencies is obtained by measuring the low 
frequency noise power after heterodyne mixing (at room temperature) the 
amplified signal from the sample against a variable frequency oscillator. 

3. Measurement of the frequency dependence of the shot noise 
in a diffusive N wire 

In this section we report measurements of the frequency dependence of the 
out-of-equilibrium (V ^ 0) noise in a normal metal wire between two N 
reservoirs (see fig. 1(a)) [3]. In such a system the current noise at finite 
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frequency is given by [12]: 



T, V) = Ak B TG (J (g(v + w)+ g(v - w)) + (1 - r,)g(w)) (6) 

The rj = 1/3 factor corresponds to the shot noise reduction due to disor- 
der (Fano factor) [12]. The shot noise corresponds to the v terms whereas 
equilibrium noise is given by v = 0. As shown by eq. (6), the total noise is 
not given by the sum of equilibrium noise (classical Johnson or quantum) 
and shot noise, even at zero frequency (w = 0). The unusual nature of this 
superposition can be emphasized by examining the fluctuations predicted 
by eq. (6) as a function of voltage for different frequencies (see fig. 4 left). 
At zero frequency (full line) , there is a transition from Johnson noise to the 
linearly rising shot noise at eV ~ fcgT ~ 2 /j,eV. At 20 GHz (dotted line), 
the fluctuations are dominated by quantum noise and do not increase from 
their value at equilibrium, until the voltage V c = Twu/e ~ 80/iV is exceeded, 
even though the condition eV > fcgT is fulfilled and the low-frequency 
fluctuations (solid line) are increasing rapidly. Only above V c does the dc 
voltage provide enough energy to increase the emitted photon noise. 

The differential noise theoretically expected for the diffusive conductor, 
under the example conditions of T = 25 mK, and a small (AV = 30^V p.p.) 
square-wave voltage modulation, are shown in the bottom left of fig. 4. The 
square- wave modulation contributes significantly to the width of the rises. 
The noise measurements are reported as a variation, AT/v, of the sample's 
noise temperature due to the voltage modulation AV. The measured values 
of AT/v for frequencies of 1.5, 5, 10, 15, and 20 GHz, taken at a mixing 
chamber temperature of T = 40 mK, are shown in fig. 4 right. While we see 
that AT/v for the low-frequency noise (circles) changes rapidly with voltage, 
approaching its linear asymptote at voltages only a few times ksT/e, the 
curves become successively broader for increasing frequency. The noise for 
the highest frequencies has a clearly different shape, displaying the expected 
plateau around V = 0. Also shown in fig. 4 right (full lines) are theoretical 
curves based on eq. (6), accounting for the finite voltage difference used, 
and for an electron temperature of 100 mK. The asymptotic value of AT/v 
has been arbitrarily scaled (since the frequency dependent system gain is 
not known to better than about 30%) to be 112 mK for each frequency, 
corresponding to the expected reduction factor of rj = 1/3 (see eq.(6)). 
Note that such a reduction of the shot noise could also be attributed to the 
heating of the electrons (hot electron regime) [13]. See ref.[3] for a detailed 
discussion of this possibility. 

Though it has not been measured, in an NS geometry the same qualita- 
tive behaviour is expected, with the scale eV replaced by 2eV. For the NS 
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Figure 4- Left: predicted dc bias voltage dependence of noise (top) for three frequen- 
cies at a bath temperature of 25 mK. The current spectral density predicted by eq. 
(6) has been converted to an equivalent noise temperature TV through the relation 
TV = S;/(4fcsG). Note that the noise is independent of bias voltage for e\V\ < Tiuj. 
The bias voltage modulation technique we employed is shown schematically and the ex- 
pected differential noise ATjv for a 30/iV p.p. modulation is also displayed (bottom). 
Right: measured differential noise for frequencies of v = lu/2tt =1.5, 5, 10, 15, and 20 
GHz, with mixing chamber temperature of 40 mK. Solid lines show the predictions of eq. 
(6) for an electron temperature of 100 mK, and accounting for the voltage modulation 
of 30 fiV p.p. 



system there might be corrections at a frequency of the order of the Thou- 
less energy, in the same way there is a signature at ~ 4E C in the voltage 
dependence of the noise measured at low frequency (see section 5). 

4. Observation of photon assited noise in an NS wire 

In this section we report measurements of low frequency noise (measured in 
the bandwidth 1.25 — 1.75GHz) emitted by an NS sample which experiences 
both a dc and an ac bias [4] . The sample is a thin Au wire between N and 
S reservoirs, as depicted in fig. 1(b). In the presence of an ac excitation at 
frequency u, the noise emitted by a diffusive wire is: 

Sj(V,T) = Ak B TGUl-i 1 ) + 2i 1 ^ 4(a)g(v + nw)\ (7) 

\ n=— oo / 

where a = 2eV ac /(huj) is a dimensionless parameter measuring the ampli- 
tude of the ac excitation voltage (note that oj denotes here the frequency of 
the ac bias; the the frequency at which the noise is measured is considered 
to be dc). This formula is valid at low (E <C E c ) and high (E » E c ) energy, 
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Figure 5. Predicted and observed shot noise of an N-S device vs. bias voltage without 
ac bias and at different powers of ac excitation at 34 GHz: (a) theory for dTpj/dV at 
T = 100 mK with no ac (dashed line) and with ac excitation at a = 1.1, 1.4, 1.7, 2.2, 2.8 
(solid lines); (b) experimentally measured cLTn /dV with no ac bias (dashed line) and with 
ac excitation powers differing by 2 dB and corresponding to the above values of a (solid 
lines); (c) rf 2 Tjv/dV 2 obtained by numerical differentiation of data in (b). (d) d 2 Tn ■ /dV 2 
vs. bias voltage at B = (solid lines) and at B — 5 T (dotted lines) with ac excitation 
at hco/2-K = 10 and 20 GHz. The curves are offset vertically by an amount proportional 
to frequency. The solid straight lines mark the expected peak locations for the N-S case 
(at B = 0): Vp ea k = hu/(2e); the dotted straight lines mark the expected peak locations 
for a normal device (at B — 5 T): V pea k = hcu/e; (e) peak location vs. frequency for 
B = and B = 5 T; the solid and the dotted straight lines are V pea k = hco/(2e) and 
Vpeak = Tiui/e, respectively. 



where G is voltage- and frequency- independent. In between, the energy de- 
pendence of the Andreev process may give rise to corrections, as is the case 
for the conductance and for the effective charge (see section 6). 

In the absence of ac excitation, the measured differential noise dT^/dV 
vs. bias voltage for the N-S device is, within 5%, twice as big as that mea- 
sured when the device is driven normal by a magnetic field of 5 T [14]. This 
is a measure of the doubling of the effective charge, which will be discussed 
more in section 5. We now turn to the noise measured in the presence of 
an ac excitation. If the transport is still elastic in the presence of ac exci- 
tation, the shot noise is expected to develop features at bias voltages such 
that qV = ±nhuj. The location of these features should be independent 
of ac power. In contrast, if the transport is inelastic, no photon-assisted 
features should occur. The derivative of the noise vs. bias voltage was mea- 
sured with ac excitation at oo/2tt = 34 GHz, at different levels of ac power. 
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Figures 5(a) and (b) show the predicted and observed derivative of the noise 
temperature vs. dc bias voltage for several levels of ac power. To see the 
features more clearly, we plot in fig. 5(c) the second derivative d 2 T^/dV 2 
obtained by numerical differentiation of the experimental data. With no ac 
excitation, d 2 T]^/dV 2 has a peak at V = 0. With ac excitation, the side- 
bands of this peak are clearly evident at V = ±nhuj/(2e). The sideband 
locations are power independent, which further argues that the structure 
is due to a photon-assisted process. The mag nitude of d 2 T N /dV 2 at V = 
displays oscillatory (roughly ~ J 2 {a)) behavior vs. ac excitation amplitude 
(not shown), which is another hallmark of a photon-assisted process. We 
note that photon assisted processes are seen clearly in SIS tunnel junc- 
tions^]. The features there are in the quasiparticle current, so the charge 
involved in that case is le. They are centered at the gap voltage V = 2A/e, 
since at low temperature pair breaking must occur for a quasiparticle to 
tunnel. 

The most convincing evidence of the photon-assisted nature of the ob- 
served effects is the dependence of the voltage location of the sideband peak 
on the frequency of the ac excitation. Measurements of the shot noise were 
made at several different frequencies of ac bias, both in zero magnetic field 
and at B = 5 T, for which the Nb reservoir is driven normal. Figure 5(d) 
shows the second derivative of the shot noise power vs. bias voltage for 
fiu)/2n = 10 and 20 GHz at B = (solid lines) and for the same device 
at B = 5 T (dotted lines), where the sample is driven normal. The solid 
and dotted straight lines are the expected peak positions for the N-S and 
normal cases, respectively. The peak locations clearly follow the theoretical 
predictions V = hto/q with q = 2e in the case of the N-S device and q = e in 
the case of the device driven normal. Figure 5(e) shows the peak locations 
for a number of different ac excitation frequencies at B = and B = 5 T. 
The solid and dotted lines are theoretical predictions with no adjustable 
parameters. 

5. Measurement of the phase dependent effective charge in an 
Andreev interferometer 

In this section we show precise measurements of the effective charge q e ff, 
and how deviations from q e ff = 2e are phase sensitive. We relate these de- 
viations to correlations of the charge transfer during Andreev reflection [5] . 
After an Andreev process, the reflected hole carries information about the 
phase of the superconducting order parameter of the S reservoir at the N-S 
interface. When two S reservoirs are connected to the same phase-coherent 
normal region, a phase gradient develops along the normal metal, result- 
ing in phase-dependent properties. In an Andreev interferometer - a device 
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Figure 6. (a) Experimentally measured effective charge q e f / for several values of mag- 
netic flux, (b) Theoretical predictions for Ec = 7.5peV and T = 43mK. The dip in q e ff 
is predicted to occur at ~ 4E C . The inset shows the theory for $ = and $ = $o/2 at 
T = 0. Note that our definition E c — hD/L 2 uses L for the full length of the normal 
region, and thus differs from the definition in Ref. [5] . 

containing a mesoscopic multiterminal normal region with a (macroscopic) 
superconducting loop, all the electronic properties are periodic with the 
magnetic flux enclosed by the loop, with a period of the flux quantum, 
$0 = h/(2e). The sample used for this experiment is depicted in fig. 1(c) 
and (d). 

From the noise measurements we deduce the effective charge, q e ff = 
(3/2)(dSj /dl); see fig. 6(a). By considering dSj/dl rather than dSj/dV 
we eliminate the trivial effect of a non-linear I(V) characteristic. At finite 
energy (E > /cgT) the effective charge reflects the charge transferred but 
also includes the effects of correlations in the transfer process. The voltage 
dependence of dSi/dl yields information on energy-dependent correlations 
between charge transfers. Figure 6(b) gives the theory results based on full 
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counting statistics. The inset shows the theory for $ = and <I> = 3>o/2 
at T = 0. The effective charge is seen in the theory to be independent of 
the phase difference at bias voltages larger than ~ 100 /iV, with signifi- 
cant phase modulation of q e ff in the bias voltage range 10 — 80 fiV. The 
maximum magnitude of the observed dip of q e ff vs. voltage is ~ 10%, and 
occurs for $ ~ $o/4- There is no dip for <E> = &o/2. For T = 0, q e ff returns 
to 2e as V — > (see inset of fig. 6). At finite temperature, q e ff goes to zero 
for eV -C ksT. This is because Si reduces to Johnson noise at V = 0. Thus, 
the decrease of q e ff at finite temperature and at very low voltages is not 
related to Andreev physics. In contrast, the dip near AE C ~ 30 fieV (with 
E c = TiD/L 2 ) is due to the energy dependence of the Andreev processes. 

The experimental results are in fairly good agreement with the theoret- 
ical predictions. As expected, there is no phase modulation of q e ff at large 
energies eV » E c , and here q e ff = 2e. At E ~ AE C , the effective charge is 
smaller for integer flux than for half-integer flux. The non-trivial energy- 
and flux dependence predicted (crossings of the different curves) is seen in 
the experiment, though the agreement is not perfect. The magnitude of the 
dip of q e ff in the data is also close to the theoretical prediction. 

To understand the origin of the dip of Sj seen for <I> = 0, we have 
also solved a generalized Boltzmann-Langevin (BL) equation. In such an 
approach correlations due to the superconductor enter only through the 
energy- and space-dependent conductivity, which gives I(V,&). Thus, the 
BL result is not complete, and we will compare its predictions to that of 
the full-counting-statistics theory, to help understand those predictions. At 
T = 0, the BL result for all flux values is simply Sf L = {2/3)2eI(V, i.e., 
q e ff = 2e at all energies. This implies that the deviation of the effective 
charge from 2e, measured and predicted by the full theory, must be due to 
fluctuation processes which are not related to single-particle scattering, on 
which the BL approach is based. We believe that the higher-order process 
which is responsible for the dip of Si is a two-pair correlation process. 
At high energies (E S> E c ) the electron-hole pair states have a length 
~ (fiD j E) 1 / 2 , shorter than L. This results in uncorrelated entry of pairs 
into the normal region. For E ~ Ec the pair size is larger, and the spatial 
overlap prevents fully random entry, suppressing Si. Suppressed shot noise 
is a signature of anti-correlated charge entry [12]. At yet lower energies (at 
T = 0) the effective charge is predicted to return to 2e; we do not yet have 
a physical interpretation of this. In any case, for the case of <E> = &o/2, the 
dip of q e ff is fully suppressed, according to the theory. This means that the 
phase gradients destroy the pair correlation effect responsible for the dip. 
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6. Theoretical approach to the effective charge 

We now turn to our theoretical approach to current noise in mesoscopic 
proximity effect structures. Our goal in this section is to explain the pre- 
dictions and the meaning of the effective charge. We shall see that the 'dip' 
in the effective charge seen for the interferometer is also seen in wires, and 
arises from similar pair correlation effects. A characteristic feature of the 
NS structures is that the phase coherent propagation of Andreev pairs in 
the normal metal is influenced by the proximity effect. One consequence is 
the so-called reentrant behaviour of the conductance of a normal diffusive 
wire in good contact to a superconducting terminal [2]. The conductance 
is enhanced at energies of the order of ~ 5E C . At higher and lower energies 
the conductance approaches its normal state value [16]. In a fork geometry, 
as discussed in the previous section, the proximity effect can be tuned by 
a phase difference between the two superconducting terminals. 

We are interested in the bias-voltage dependence of the current noise in 
a diffusive NS structure. There is, however, a simple energy dependence re- 
sulting from the proximity-induced energy-dependent conductivity. In con- 
trast, the correlations of interest are not due to this energy-dependence of 
the conductivity. They can be distinguished in the following way. We note 
that for the average transport properties, i.e. the differential conductance, 
the kinetic equation takes the very simple form [16, 17, 18, 19] 



The energy- and space-dependent conductivity cr(x, E) includes the prox- 
imity effect, and /t(x, E) = 1 — /(x, E) — /(x, — E) is the symmetrized 
distribution function. Due to the induced superconducting correlations, a is 
enhanced above its normal state value ajy and its energy- and space depen- 
dence is obtained from the spectral part of the Usadel equation [16, 20]. For 
the geometry in fig. 7 the spectral conductance is G~ 1 (E) = J dx/a(x, E). 
Note that G(E) ^ even for E <C A as a consequence of the proximity 
effect. The current for a given bias voltage V and temperature T is then 
given by 



where fjf(E,V,T) is the symmetrized distribution function in the normal 
metal reservoir and we have accounted for the boundary condition fj, (E <C 
A, V, T) = at the superconducting terminal. 

The form of the kinetic equation (8) suggests that electrons and holes 
(i.e. positive and negative energy quasiparticles) obey independent diffusion 
equations, which are only coupled through the boundary condition fj, = 



W(x,£)V/ T (x,£) =0. 



(8) 



1 f + OO 

I(V,T) = - / G(E)fJ*(E,V,T)dE 
2e J-oc 



(9) 
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at the superconducting terminal. Thus, we may try to apply the semiclas- 
sical Boltzmann-Langevin (BL) approach [21, 22]. The only modification is 
that we have to account for an energy- and space-dependent conductivity. 
For that purpose it is convenient to introduce the characteristic potential, 
defined as the solution of the equation : 

X7a(x,E)Vis(x,E) = (10) 

with the boundary condition that v = 1 at the normal terminal and v = at 
the superconducting terminal. The solution for our quasi-one dimensional 
geometry is : 

The current noise can then be expressed in the familiar form : 

Sf L {V) = A J dE J dxa(x,E) (Vu(x,E)) 2 f(x,E)(l - f(x,E)) (12) 

where the distribution function is given by f(x,E) = v(x, E)f N (E,V,T). 
As a result we find for the noise at zero temperature (i.e., fjf(E, V, T = 0) = 
-sign(eF) for \E\ < \eV\ and zero otherwise) : 

Sf L (V) = ^I(V) (13) 

where the current is given by /(V) = (l/2e) J^ v dEG(E). Thus, we find 
that the current noise depends in a non linear fashion on the voltage. The 
nonlinearity is given by the I{V) characteristic. However, this dependence is 
in some sense trivial, since the only way the electron-hole coherence enters 
is through the energy dependent conductivity. 

We note that the doubling of the shot noise in comparison to the nor- 
mal case results from the energy integration from —eV to +eV, instead of 
the interval between (i.e., the Fermi energy) and eV, as it would be in 
the normal case. On the other hand, the average 'noisiness' (coming from 
the spatial integral in eq. (12) alone) is the same in the normal and the 
superconducting cases [22]. Thus, the doubling of the noise needs not be 
interpreted as a direct consequence of the doubled charge transfer involved 
in an Andreev reflection process. It reflects the particle-hole symmetry in 
the superconducting terminal. Nevertheless we adopt below the notion of 
an effective charge, since it is a convenient measure of the deviation from 
the independent electron fluctuations. 

These arguments can also be used to explain the experiments on photon- 
assisted noise described previously. In the presence of an ac voltage of fre- 
quency to the electron distribution in the normal terminal acquires side- 
bands, i.e., additional steps at energies ±nhuj. The noise of the diffusive 
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Figure 7. Transport characteristics of a proximity wire. The inset shows the layout. 
The relevant energy scale is the Thouless energy E c = hD/L 2 . In the main plot the 
differential conductance, the differential noise dSi(V) / dV , and the effective charge as 
defined in (14) at zero temperature are shown. 



wire depends essentially on a superposition of left and right distribution 
functions differing by the voltage eV in the normal state and by 2eV in the 
NS case. It is clear that the noise properties change as a function of voltage, 
when sideband features match the other steps in the distribution function. 
As a consequence, photon- assisted steps occur when the voltage matches 
nfitu/e in the normal case and nfiLojle in the superconducting case. This is 
what is observed in the experiments (see section 4). 

A correct calculation of the noise requires that we go beyond the in- 
dependent electron- and hole-fluctuations in the Boltzmann-Langevin ap- 
proach. This can be accessed by the extended Green's function approach 
[23]. To describe the fluctuations, we utilize the results of the full counting 
statistics and define an effective charge 

3 as f (y,r) 

Qeff\V,T) — 2~df(y^T) ' (14) 

which takes the value 2e for the Boltzmann-Langevin result (see eq. (13)). It 
follows that the energy dependence of the effective charge gives information 
about the correlated electron-hole fluctuation processes. 

In order to illustrate these considerations, we show in fig. 7 results for 
the conductance dl/dV, the differential noise dSi(V)/dV, and the effective 
charge at zero temperature for a one-dimensional diffusive wire between 
a normal and a superconducting reservoir. Note: q e ff is proportionnal to 
dSi/dl = (dSj/dV)(dI/dV)~ 1 . These results were obtained by a numeri- 
cal solution of the quantum-kinetic equation. The differential conductance 
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shows the well-known reentrance (peak) behaviour [16]. At low and high 
energy the conductance approaches the normal state value Gat. At inter- 
mediate energies of the order of the Thouless energy the conductance is en- 
hanced by ~ 15% above Gn- A similar result is found for the interferometer 
when $ = [5] . We observe that the differential noise dS/ dV has a roughly 
similar energy dependence [23], although it is quantitatively different. The 
deviation of the voltage-dependent effective charge from 2e demonstrates 
the energy dependence of the higher order correlations, which are not con- 
tained in the independent electron-hole picture, the BL approach. Around 
E ~ AE C the effective charge is suppressed below 2e, showing that the 
higher order correlations result in a reduced noise in comparison to the BL 
case of uncorrelated electron and hole fluctuations, for which q e ff = 2e. As 
V — > 0, q e f f approaches 2e again. 

The physics of this effective charge is clarified if we consider the full 
counting statistics [23, 24, 25, 26], instead of the current noise only. In full 
counting statistics, we obtains the distribution of tranfered charges, which 
clearly contains direct information on both the statistics and the nature 
of the charge carriers. There it follows that all charge transfers at subgap- 
energies occur in units of 2e [27]. However, this does not necessarily result 
in a doubled effective charge using our definition of q e ff- The effective 
charge also includes the effect of correlations between the different charge 
transfers. Our work in section 5 shows that these correlations are phase- 
dependent. Nevertheless, in the limit of uncorrelated transfer of Andreev 
pairs, the effective charge at E S> ksT is simply 2e. 

7. Conclusion 

In this paper we have discussed how Andreev Reflection affects shot noise 
of mesoscopic NS structures. High frequency measurements provide a very 
powerfull tool since they are very sensitive, as is essential to the measure- 
ment of the small voltage- and flux dependence of the effective charge, 
and these measurements access an energy domain in which interesting phe- 
nomena occur, when hu > eV,ksT. These techniques are very promising 
for investigation of even more subtle quantities like cross-correlations in the 
noise or higher moments of the current fluctuations [28] . The measurements 
have revealed the existence of correlations in pair charge transfers which 
are not accessible through conductance measurements. The full counting 
statistics method gives access to the full distribution of the charge trans- 
fers, and is essential for understanding the physics of such NS structures. 
This method sheds light on the correlations revealed by the shot noise, and 
allows as well the investigation of cross-correlations and higher moments. 
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